Astronomy & Astrophysics manuscript no. 1518-astroph 


© ESO 2008 


August 17, 2008 





Multi-telescope timing of PSR J1 51 8+4904 

G. H. Janssen^ B. W. Stappers^'^'^, M. Kramer^, D. J. Nice^, A. Jessner^, I. Cognard^, and M. B. Purver^ 

' Astronomical Institute "Anton Pannekoek", University of Amsterdam, Kruislaan 403, 1098 SJ Amsterdam, The Netherlands; 

e-mail: geitiina@science.uva.nl 
^ Stichting ASTRON, Postbus 2, 7990 AA Dwingeloo, The Netherlands; e-mail: Ben . StappersOmanchester .ac.uk 
^ University of Manchester, Jodrell Bank Observatory, Macclesfield Cheshire, SKI 1 9DL, UK 

Physics Department, Bryn Mawr College, Bryn Mawr, PA 19010, USA 
^ Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53 121 Bonn, Germany 

* Laboratoire de Physique et Chimie de I'Environnement, CNRS, 3A Avenue de la Recherche Scientifique, F-45071 Orleans, Cedex 2, France 
Received/ Accepted 

ABSTRACT 

Context. PSR J15 18+4904 is one of only 9 known double neutron star systems. These systems are highly valuable for measuring the masses of 
neutron stars, measuring the effects of gravity, and testing gravitational theories. 

Aims. We determine an improved timing solution for a mildly relativistic double neutron star system, combining data from multiple telescopes. 
We set better constraints on relativistic parameters and the separate masses of the system, and discuss the evolution of PSR J15 18+4904 in the 
context of other double neutron star systems. 

Methods. PSR J1518+4904has been regularly observed for more than 10 years by the European Pulsar Timing Array (EPTA) network using 
the Westerbork, Jodrell Bank, Effelsberg and Nanjay radio telescopes. The data were analysed using the updated timing software tempo2. 
Results. We have improved the timing solution for this double neutron star system. The periastron advance has been refined and a significant 
detection of proper motion is presented. It is not likely that more post-Keplerian parameters, with which the individual neutron star masses and 
the inclination angle of the system can be determined separately, can be measured in the near future. 

Conclusions. Using a combination of the high-quality data sets present in the EPTA collaboration, extended with the original GBT data, we 
have constrained the masses in the system to njp < 1.17 Mq and m^. > 1.55 Mq (95.4% confidence), and the inclination angle of the orbit to 
be less than 47 degrees (99%). From this we derive that the pulsar in this system possibly has one of the lowest neutron star masses measured 
to date. From evolutionary considerations it seems likely that the companion star, despite its high mass, was formed in an electron-capture 
supernova. 

Key words, stars: neutron - pulsars: general - pulsars: individual (PSR J15 18+4904) 



1. Introduction 



PSR J15 18+4904 was discovered in the Green Bank Northern 
Sky Survey ( |Sayer et al.|1997 1. It is one of only 9 double neu- 
tron star systems (DNSs) known. The pulsar orbits its com- 
panion neutron star in 8.63 days, and its 40 ms spin period, 
together with the derived low surface magnetic field, is typical 
of a mildly recycled pulsar. |Nice et aL] ([1996 ) first described 
the system and have akeady shown that the space velocity is 
probably quite low. They also measured the periastron advance 



from which the total system mass was estimated. Thorsett & 



Chakrabarty (1999i, using a Bayesian analysis with no con- 
straints on the orientation of the orbit, found the masses of 
the pulsar and its companion to be Wp = 1.56^j 2q Moand 
OTc = 1-05^Q [! Mq (95% confidence). An improved timing so- 



lution was presented by Hobbs et al. ( 2004j ), although not dis- 
cussed in detail. 



DNSs are presently the best available tool for testing 
strong-field gravity eff'ects. The PSR J1518+4904 orbit is only 
mildly relativistic, making it of Umited use for tests of gravita- 
tional theories. However, any constraints on its post-Keplerian 
(PK) parameters, the inclination of its orbit, or the masses of 
the pulsar and its companion are highly valuable for studies of 
the evolution of these systems (e.g. Podsiadlowski et al.|2004 
Kalogera et al.|2007| l. 

Like PSR J1811-1736 (|Corongiu eTlIl [20071) 



PSR J15 18+4904 has a rather wide orbit. This could indicate 
that the evolution into a DNS system has been slightly different 
than most other systems, which generally have tight orbits with 
periods of several hours. However, PSR J1518+4904 does 
follow the recently discovered spin period-eccentricity relation 



(Faulkner et al. 2005 Dewi et al. 2005 i, suggesting that the 
evolution cannot be too different from the tight DNSs. 

PSR J 15 18+4904 has been observed regularly using the 
four 100 m class radio telescopes in Europe. This has allowed 
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us to select high-quality data resulting in the best possible anal- 
ysis of the system parameters to date. For completeness, we 
have also included data from telescopes at Green Bank as pre- 



sented in Nice et al. ( 1996 1999j ). We present the new timing 
solution and discuss the limits our solution puts on the incli- 
nation and the masses of the neutron stars in the system. We 
discuss the future prospects of detecting multiple PK parame- 
ters, and we related our results to the evolution of DNSs. 

2. Observations and data analysis 

This research is part of the EPTA network. The extensive 
collection of observations from Westerbork, Jodrell Bank, 
Effelsberg and Nanfay observatories has enabled us to select 
the best combinations of data sets with the best characteristics 
for various purposes. For example, high quality TOAs and long 
time spans permit better measurement of astrometric param- 
eters, while using multiple frequencies allows us to precisely 
measure the dispersion measure (DM) and monitor possible 
DM variations. Individual properties of the data sets are pre- 
sented in Table[T[ and a description of the observing systems is 
given below. 

2.1. Westerbork 

PSR J 15 18 -1-4904 has been observed approximately monthly 
since 1999 at the Westerbork Synthesis Radio Telescope 
(WSRT) with the Pulsar Machine (PuMa; [\bfl"te et al.|2002| l, at 
frequencies centred at 840 and 1380 MHz, and since 2007 also 
occasionally at 2300 MHz. The sampling time for most obser- 
vations was 102.4 jUS and the bandwidth used was 8x10 MHz 
where each 10 MHz band was split into 64 channels. The 1380 
MHz data taken after September 2006 used 80 MHz of band- 
width spread in 8 steps of 10 MHz over a range of 160 MHz. 
The data were dedispersed and folded off-line, and then inte- 
grated over frequency and time over the whole observation du- 
ration to get a single profile for each observation. Each profile 
was cross-correlated with a standard profile (Fig. [T}, obtained 
from the summation of high signal-to-noise (S/N) profiles, to 
calculate a time of arrival (TOA) for each observation. These 
were referred to local time using time stamps from a H-maser 
at WSRT. The TOAs were converted to UTC using global po- 
sitioning system (GPS) maser offset values measured at the ob- 
servatory. For this research, we only used observations with 
TOA errors < 15 ^s, with an average around 8 /is. 

2.2. Jodrell Bank 

At Jodrell Bank, PSR J1518+4904 was observed since its 
discovery approximately every two weeks using the 76-m 
Lovell telescope. Observations were made at centre frequen- 
cies around 410, 610 and 1400 MHz. All receivers were cryo- 
genically cooled, providing left-hand and right-hand circularly 
polarized signals. These signals were fed into an analogue fil- 
terbank system where the polarisations were detected, filtered, 
digitised at appropriate sampling intervals and incoherently 
dedispersed in hardware. The sampling time was chosen to 
match the dispersion smearing across one filterbank channel. 




Fig. 1. Standard profiles for the three different frequencies used 
at Westerbork. From top to bottom: 840 MHz, 1380 MHz, 2300 
MHz. The effective time resolution (102.4 jUs) is similar to the 
size of one bin (160 jUs). 



At 410 MHz, a 2 X 32 X 0.0312 MHz system was used, at 610 
MHz we employed a 2 x 6 x 0.1250 MHz filterbank, while the 
vast majority of our data were recorded with a 2 x 32 x 1 MHz 
system at 1400 MHz. 

The resulting dedispersed timeseries were folded on-line 
with the topocentric pulsar period and finally written to disc. 
In the ofi^-line reduction, the two polarisations were summed 
to form total-intensity profiles. A standard pulse template was 
used at each frequency to determine the TOA. During this pro- 
cess, TOAs were referred to the local H-maser time-standard 
and already corrected to UTC using information obtained via 
the GPS. In our analysis, we rejected TOAs with estimated un- 
certainties exceeding 40 /is (or 60 /is for data prior to Jan 2000, 
respectively). With an average TOA error for the final data set 
of 18 /vs, the precision of this incoherently dedispersed data set 
is modest, but it provides the longest time baseline in our anal- 
ysis. The Jodrell Bank data were therefore used as a reference 
to calculate observatory-related off'sets between the other data 
sets used. 



2.3. Effelsberg 

Observations in Effelsberg were made at least once a month at 
1410 MHz and at irregular intervals at 850 MHz. At 1410 MHz 
a tunable HEMT receiver was used achieving typical system 
temperatures of 30-40 K depending on weather conditions and 
telescope elevation. 

At the highest frequency a total bandwith of 55 MHz was 
available from the coherently dedispersing backend. For each 
polarisation the band was split into four subbands which them- 
selves were subdivided into eight digitally sampled channels. 
Each of these 32 bands per polarisation was coherently dedis- 
persed by programmable digital filters and then recombined 



with the approriate channel dispersion delay (Backer et al. 



1997 1 and synchronously summed up with the current topocen- 
tric period. 
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Table 1. Individual data sets. 





Effelsberg* 


Jodrell Bank 


Nan§ay* 


Westerbork 


Green Bank* 


All combined 


Number of TOAs 


71 


292 


145 


126 


382 


1016 


Time span (MJD) 


52481 -54166 


49797 - 53925 


53307 - 54200 


51389-54337 


49670 - 52895 


49670 - 54337 


R.M.S. individual data set" (jUs) 


3.3 


10.7 


3.3 


5.2 


8.2 


6.0 


Observed Frequencies (MHz) 


1400 


400, 600, 1400 


1400 


840, 1380, 2300 


350, 370, 575, 800 





The best timing solution acquired from all data sets is presented in Table |2] (a) R.M.S. given require reduced = 1 for each individual data 
set. (b) R.M.S. is the result of a fit with proper motion (Green Bank, Effelsberg, Nan^ay) and dispersion measure (Effelsberg, Nanjay) fixed. 



A TOA was calculated for each average profile obtained 
during a 5-10 minute observation. During this process, the ob- 
served time-stamped profile was compared to a synthetic tem- 
plate, which was constructed out of 3 Gaussian components 
fitted to a high S/N standard profile ( [Kramer et al. 1 1 998| [T999 1 . 
This template matching was done by a least-squares fitting of 



the Fourier-transformed data (Taylor 19921. Using the mea 



sured time delay between the actual profile and the template, 
the accurate time stamp of the data provided by a local H-maser 
and corrected off-line to UTC using recorded information from 
GPS satellites, the final TOA was obtained. 

Originally, multiple TOAs were calculated from one ob- 
servation, but better results were acquired from calculating one 
TOA per observation. We have used data obtained since August 
2002, having TOA errors smaller than 20 /is, on average 5yus. 



2.4. Nangay 

The pulsar has been observed weekly in Nanfay since late 
2004. Equivalent to a 93-m dish, the Nangay Radio Telescope 
and the BON (Berkeley-Orleans-Nangay) coherent dedispersor 
were used for typical integration times of 45 min. Coherent 
dedispersion of a 64 MHz band centred on 1398 MHz is car- 
ried out on sixteen 4 MHz channels using a PC-cluster. The 
Nangay data are recorded on a UTC(GPS) time scale built at 
the analogue to digital converter by a Thunderbolt receiver 
(Trimble Inc.). Differences between UTC and UTC(GPS) are 
less than 10 ns at maximum, and therefore no laboratory 
clock corrections are needed. One TOA was calculated from 
a cross-correlation with a pulse template for each observation 
of ~ 45 min. Except for 24 low S/N profiles that were ex- 
cluded from the used sample, all TOAs have errors less than 
15 /is, but generally around 5 /is. Since the data span of Nangay 
Observatory covers only a few years, the proper motion of 
PSR J 15 18-1-4904 was fixed in determining the best timing so- 
lution for this data set alone. 



2.5. Green Bank 

PSR J 15 18-1-4904 was observed at Green Bank using the 140 
Foot (43 m) Telescope and the 100 m Green Bank Telescope 
(GBT). Observations with the 140 Foot Telescope were made 
every few weeks from November 1994 to July 1999 at 575 and 
800 MHz. Intensive campaigns, giving full coverage of the or- 
bit over several days, were made with the 140 Foot Telescope at 
370 MHz in March/April 1995, August/September 1995, May 



1998, and December 1998. An intensive campaign was made 
using the GBT at 370 MHz in August/September 2003. 

The same data acquisition system and nearly identical ob- 
servational procedures were employed at the two telescopes. 
Details are given in |Nice et al.|p996| . Briefly, data were col- 
lected by the Spectral Processor, a Fourier transform spectrom- 
eter, which synthesized 512 spectral channels across a pass- 
band of 40 MHz and folded the power measurements in each 
channel at the topocentric pulse period to produce pulse pro- 
files. Observation times ranged from a few minutes to sev- 
eral hours on any given day. During the course of an observa- 
tion, data were integrated over intervals of 2 minutes (140 Foot 
Telescope) or 0.5 minutes (GBT). Data were separately col- 
lected in two polarizations. Off-line, each integration was pro- 
cessed by summing polarizations, dedispersing and summing 
spectral channels to produce a single pulse profile, and deriv- 
ing a TOA using conventional procedures. Observations of up 
to one hour were averaged together by computing TOAs for in- 
dividual integrations and averaging them into a single effective 
TOA for the hour-long span. Time stamps were provided by a 
local H-maser and corrected off-line to UTC using recorded in- 
formation from GPS satellites. Arbitrary offsets were allowed 
in the timing solution between sets of TOAs collected at differ- 
ent frequencies, and between two sets of GBT TOAs collected 
with different receiver polarization settings. 

2.6. Combination of data sets 

Combining data sets obtained by different observatories will 
have many advantages. Using a longer time span sampled with 
more data is usually better for timing in general. The four tele- 
scopes now in use for EPTA timing purposes all contribute 
differently to the total picture of high-precision timing. The 
Westerbork, Effelsberg and Nangay data sets are of very high 
quality, and adding the Jodrell Bank long-term data which has 
only slightly lower average R.M.S. provides us with both a 
good coverage in time as well as in orbital phase for this pulsar. 
The addition of the original GBT data further complements our 
set of measurements of TOAs from the system. 

When combining data from multiple telescopes, all with 
different observing and operating systems, it is important to 
account for all extra time coiTections needed, apait from the 
usual coiTection from arrival times at individual telescopes to 
aiTival times in TAI at the solar system barycentr^ (SSB). 
Observatory-related time delays that are not accounted for (i.e. 
unmodelled cable delays) as well as a different approach in cal- 

' ftp://ssd.jpl.nasa.gov/pub/eph/export/DE405/de405iom.ps 
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Table 2. Timing solution using data from Jodrell Bank, 
Westerbork, Effelsberg, Nangay and Green Bank. 

Fit and data-set 



Pulsar name 

MJD range 

Number of TO As 

Epoch (MJD) 

R.M.S. timing residual (jis) . 

Weighted fit 

Reduced value 



J15 18+4904 

49670.7—54337.7 

1016 

52000 

6.05 

Y 

1.117 



Measured Quantities 



Right ascension, a (J2000) 

Declination, S (J2000) 

Pulse frequency, v (s"') 

First derivative of pulse frequency, v (s^^) 

Dispersion measure, DM (cm^^pc) 

Proper motion in RA, yu„ (masyr"') 

Proper motion in DEC, yUj (mas yr"') .... 

Orbital period, Pb (d) 

Epoch of periastron. To (MJD) 

Projected semi-major axis of orbit, x (It-s) 

Longitude of periastron, oiq (deg) 

Orbital eccentricity, e 

First derivative of orbital period, Pb 

First derivative of x, x 

Periastron advance, cj (degyr"') 



15''18"lff.799084(16) 

-^49°04'34'.'251 19(16) 

24.4289793809236(3) 

-1.62263(12)xl0-" 

11.61139(8) 

-0.67(4) 

-8.53(4) 

8.6340050964(11) 

52857.71084163(17) 

20.0440029(4) 

342.554394(7) 

0.24948451(3) 

2.4(22)xl0-" 

-l.l(3)xlO-"' 

0.0113725(19) 



Derived Quantities 



logio(Characteristic age, yr) 10.38 

log,o(Surface magnetic field strength, G) 9.03 

Mass function (Mq ) 0.1 15988 

Total mass, Mt (Mq ) 2.7183(7) 

Total proper motion, fij (mas yr"') 8.55(4) 

DM distance (pc) 625!* 



Assumptions 


Clock correction procedure 


TT(TAI) 


Solar system ephemeris model 


DE405 


Binary model 


DD 


Model version number 


5.00 



Figures in parentheses are the nominal Icr tempo2 uncertainties in the 
least-significant digits quoted. All time-dependent variables refer to 
the same epoch. The DM distance is estimated from the jCordes &| 
|Lazio| ( (2002l > model. 



culating TOAs (template differences) will result in a time offset 
between sets of residuals from the different telescopes. The new 
timing software package tempo2 ( Hobbs et al.|20 06^ is capable 
of fitting for constant time offsets, and these so-called "jumps", 
which are on the order of 0. 1 ms, have been used as extra pa- 
rameters in the timing model. A more detailed study on this 
subject will be presented in a forthcoming paper. 



3. Results 

The best overall timing solution for PSR J1518+4904, ob- 
tained by combining data from Jodrell Bank, Westerbork, 
Effelsberg, Nangay and Green Bank, is presented in Fig. |2] and 
Tables [T] and |2] Before combining the data sets from the dif- 



ferent observatories, the errors on the TOAs were scaled with a 
constant factor to have each individual data set return a reduced 

Having multiple sets of TOAs allows us to select different 
subsets of TOAs for refining different parameters. When mea- 
suring the astrometric and rotational parameters like position, 
proper motion and spin frequencies, it is best to use a combina- 
tion covering a data span as long as possible, in our case using 
all available data sets. Even though the standard deviation of 
the residuals from zero was, on average, somewhat larger, this 
yields the best possible determination of these parameters. 

For short-period binaries, it can sometimes be better to se- 
lect only the data set(s) with the highest precision TOAs to de- 
termine the binary parameters more accurately. We have tried 
several combinations of the best TOAs, i.e. using Nan9ay and 
Westerbork, or Nangay, Westerbork and EflFelsberg. However, 
we found the best parameter measurements are derived from 
timing solutions which combine all data sets, not just those 
with the most precise TOAs. The total data set has broader or- 
bital coverage compared to any combination of subsets having 
the most precise TOAs, and also covers a longer time span, im- 
proving the measurement precision of parameters that undergo 
secular changes. The timing parameters derived from various 
combinations are all in agreement with each other, showing 
that there are no large systematic effects in the separate data 
sets. 



3.1. Dispersion measure variations 

Although the majority of all observations were done at frequen- 
cies around 1400 MHz, Westerbork, Jodrell Bank and Green 
Bank have also observed PSR J1518+4904at other frequen- 
cies, see Table [T] allowing us, in principle, to detect DM vari- 
ations. We applied the "stridefit" plugin of tempo2 which fits 
small segments of time for DM while keeping all other param- 
eters fixed. Because the observations at additional frequencies 
are not distributed evenly across the full time span, it was not 
possible to apply this method continuously. However, there is 
no evidence of large variations or trends present in the disper- 
sion measure. 

In a similar way as explained above, we have tried to im- 
prove our timing solution by fixing the best DM found from 
the solution using all available TOAs, and fitting again using 
only the higher precision 1400 MHz data from Jodrell Bank, 
Westerbork, Nangay and Effelsberg. Again, having the data set 
as extended as possible resulted in the better determination of 
the parameters. 

3.2. Parallax 

Even though the DM is quite small, the error in the DM-derived 
distance is still dominating the significance of any velocity or 
orientiation measurement. For a relatively nearby system like 
PSR J1518+4904 , it is worthwhile investigating the possibility 
of detecting a parallax signature in the timing residuals to be 
able to acquire an independent, and perhaps even better, dis- 
tance estimate. 
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The expected amplitude in the timing 



to parallax (from Lorimer & Kramer 2005 

f cos/3/2cd 



residuals due 
§8.2.5) is 



0.9jjs, where / is the distance between the Earth 
and the SSB, /3 the ecliptic latitude of the pulsar, c the speed of 
light and d the distance to the pulsar. Unless the TOA measure- 
ments improve significantly, we cannot expect to measure the 
parallax for at least another 600 years. 

3.3. Proper motion 

As part of the astrometric fit, we have made a significant de- 
tection of the proper motion of this system, in right ascension 
yUfj, = -0.67(4), and in declination /jg = -8.53(4) masyr"'. The 
total proper motion is fij - 8.55(7) masyr"' which, using the 
DM-derivec|^distance of 625^g3 pc, corresponds to a transverse 
velocity of just 25(4) km s"' . 

3.4. Sliklovsl<ii effect 

The transverse motion of the system results in an increasing 
projected distance of the pulsar to the solar system barycen- 
tre, which affects any observed change in periodicities in the 
system ( |Shklovskii|1970| l. 

For this pulsar, with a period of 40.93 ms, a total proper 
motion measured to be 8.55 mas yr"' and distance estimated 
at 625 pc, the Shklovskii term in the period derivative will be 
as large as 20% of the observed P. However, as will be shown 
in ^3.6.3 this effect can be cancelled out by the acceleration 
towards the Galactic disc. The observed period derivative of 
the pulsar is therefore likely to be close to the intrinsic value. 
The estimate of the characteristic age and the magnetic field 
will not deviate significantly from the tempo2 calculated values 
as presented in Table |2] 

3.5. Secular perturbation of x and d» 

Because of the significant proper motion measurement, it is in- 
teresting to investigate any geometric effects that are likely to 
result from changes in the line of sight towards the system. 
For example, the varying projection of the orbit as the sys- 
tem moves through space can result in an apparent change in 
the semimajor axis of the orbit, i, and a perturbation of the 



observed periastron advance, 6d) . Following Kopeikin ( 1996 1, 
but writing in terms of position angles, these can be expressed 
as: 



X — X jdj cot / cos(0^, - D.), 
6cj - -fij CSC / sin(0^ - Q), 



(1) 
(2) 



where and Q. represent the direction of proper motion and 
the position angle of the ascending node respectively (mea- 
sured North through East). When it can be shown that the 
proper motion is the dominant effect on the observed i , Eq. [T] 
can be used to derive an upper limit on the orbital inclination: 



tan / < fij 



^ The DM distance is estimated from the 
model. 



Cordes & Lazio 



(3) 
j2002| 



where we have used | cos(0^ - Q.)\ < 1 . In this system, the 
effect of gravitational wave (GW) emission on the change of 
semi-major axis will be negligible, as it is expected to be only a 
fraction of the GW-effect on , which is already undetectable 



(see 5 3.6.3 1, so that we can apply the above formula. 

We were able to significantly measure a change in the pro- 
jected semi-major axis: i = -1.1(3) x 10"'^. Taken at face 
value, our measurement corresponds to an upper limit on the 
inclination of 69 degrees, assuming 0° < / < 90°. A com- 
plete analysis of the interpretation of i, including all orbital 
orientations and including the effects of its covariance with the 



relativistic parameter y, will be presented in 5 4. 1 below. 



3.6. Post-Keplerian parameters 

As mentioned above, the wide orbit of PSR J 15 18-1-4904 
means that it is unlikely that another PK parameter, besides d> , 
will be detected in the near future. However, considering that 
Pb is the PK parameter with the strongest dependence on the 
total observing time span T and Pb ( [Damour & Taylor||1992[ 
Table II), it will probably be measured sooner than any of the 
other PK parameters. Placing an exact value on the timescale it 
will take to measure any further PK parameters is difficult, due 
to the covariances with other parameters and the possibility that 
kinematic terms may come to dominate the timing solution. In 
general, measuring any further PK parameters in this system 
will be challenging and will require observation durations that 
at least double the existing observation span of 12.8 years. 

We consider each of the PK parameters below. For this sys- 
tem we neglect any contributions to the PK parameters due to 
effects other than general relativity (GR) or proper motion. As 
spin-orbit coupling and tidal effects are only expected to be of 
influence for a main-sequence companion, and any significant 
third body in the system would have been detected in the timing 
residuals, we consider this a valid assumption. 



3.6.1 . Advance of periastron: to 

We measure a value of w= 0.0113725(19) degyr^^ for the 
periastron advance of the orbit. This agrees with, and im- 
proves upon the alre ady published values of 0. 01 1 1(2) deg yr"' 
0.0113(1) degyr- 



1996 19991 and 0.01138(4) 



( Nice et al. 
degyr-i ( |Hobbs et al.|200 4'). 

Our measured cii , combined with the mass function result- 
ing from the orbital parameters, gives a very precise total mass 
of the system of 2.7183(7) Mq . This value is calculated assum- 
ing GR is the correct theory of gravity and the only measurable 
cause of periastron advance. As can be seen from Eqj2] and 
shown in more detail in Fig.|7] this is likely justified, as we find 
the contribution to d> from the proper motion for the possible 
geometric configurations of the system to be of the same order 
as the measurement uncertainty. 



3.6.2. Sliapiro delay parameters r and s 

For suitable inclination angles, Shapiro delay would be readily 
detected in this system. However, we have not seen the signa- 
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ture of Shapiro delay in the data, which indicates that the orbit 
has a low inclination. We use this to set a limit on the incli- 
nation of the orbit by generating a map of space resulting 
from fitting the mass of the companion and inclination angle in 
various combination^ The resulting probability contours are 
shown in Fig. [3] along with the constraint on the total mass 
from our (1) measurement. In our case the contour plot results 
in upper limits only and there is a 99% chance that sin / < 0.73, 
corresponding to an upper limit on the inclination angle of 47 
degrees. Combining this limit with the d) measurement and the 
mass function gives a lower limit on the companion mass of 
1.29 Mq. 

3.6.3. Decay of the orbit: 

Another effect of the wide orbit of the system is that the de- 
cay of the orbit due to gravitational wave damping will be 
very small. The relative motion effects, due to the proper mo- 
tion and the acceleration in the Galactic potential, wil l there- 
fore completely dominate observed changes in P^, (e.g. Bell & 
Bailes 1996[ Nice & Taylor 19951. Even with the relatively low 
proper motion measured for this system, the Shklovskii effect 
( |Shklovskii|1970| l on the orbit will be 



c d 



9.5 X 10 



-14 



(4) 



where we have used the aforementioned distance of 625 pc. 

The acceleration toward the disc, though of opposite sign, 
has roughly the same magnitude. The net bias of P\, is probably 

in the range APb (1 - 3) x 10 depending on the distance 

to the pulsar The uncertainty in this bias is much larger than 
the expected P^, due to general relativity, /"b.cR 1.2 x 10"'^. 

Our timing solution yields an upper limit of < 2.4 x 
10 This is considerably larger than that expected from rel- 
ativistic decay of the orbit or from kinematic effects in the 
Galaxy. Although P^ may be the first other PK parameter to 
be measured as the data set is extended, the combination of the 
three effects prevent the use as either an independent distance 
estimation or the separation of the masses of the system. 

4. Discussion 

4.1. Masses 

For any given theory of gravity, measuring PK parameters will 
put constraints on the masses and the inclination of the system, 
depending only on the individual masses and the Keplerian 
parameters (e.g. Stairs 2003 | l. Figure |4] shows all the current 
restrictions we were able to derive from our timing solution 
(Table |2]i. The improved, highly significant (1> measurement 
gives us a very accurate determination of the total mass. From 
the non-detection of Shapiro delay parameters it is clear that 
the inclination of the orbit is quite low. Furthermore, the i mea- 
surement presented in §3.3, taken at face value, confirms that 
the system is at low inclination angles. 



The relativistic time-dilation/gravitational redshift parame- 
ter, y, is only measurable when the periastron advance is large 
enough to decouple its effect on the timing residuals from 
the measurement of the semimajor axis and its time deriva- 



tive (g landford & Teukolsky|1975| [Manchester & Taylor] 1977 
|Damo ur & Tay lor||1992| l. Changing o) by only a few degrees 
will take several hundred years, and a detection of y can there- 
fore not be expected in the near future. Moreover, i turns out 
to be highly covariant with y in the timing analysis. This makes 
an independent measurement of y impossible and obscures the 
interpretation of the apparent measurement of i . 

To clarify the influence of all the PK perturbations on our 
timing data, and to calculate accurate, refined values of the stel- 
lar masses, we performed a comprehensive, self-consistent tim- 
ing analysis which simultaneously incorporated all relativistic 
and kinematic phenomena in the timing solution. 

Our primary goal was to measure or constrain the values 
of Wp and mc. This procedure also finds constraints on the two 
angles which describe the orientation of the orbit, / and O. For 
convenience, we always represent the latter as an offset from 
the proper motion position angle, i.e. 0^, - Q. This angle runs 
between -180° and 180°, while inclination runs between 0° 
and 180°. 

The four degrees of freedom in this problem (nip, trie, U and 
0^ - O) are reduced to three by noting that the two masses 
and the inclination are related via the Keplerian mass function 
equation: 



(mc sin 0"' 
(Mt)2 " 7^ 



(5) 



^ This analysis used the tempo2 plugin m2sini. For explanation of 
this plugin, see Ord et al. ( 2006| . 



where all quantities on the right side of the equation are mea- 
sured to high precision. We analysed a grid of timing solutions 
in a three-dimensional parameter space following the approach 
as explained in detail in |Splaver et al.| ( |2002| [2005 | l. For the 
three variables, we used Mj, cosi, and 0^ - Q, and we as- 
sumed a uniform prior in each of these. For randomly oriented 
binary systems, cos / and 0^ - Q each follow uniform distri- 
butions, and since cii is well known, only a small range of Mt 
values need be considered, so a uniform prior is acceptable for 
this variable as well. For each point in our three dimensional 
grid, we used Eq. |5]to calculate the value of and the corre- 
sponding mp. We then used the masses and orientation angles 
to calculate the kinematic perturbation parameters (Eqs. [T| [2]i 
and relativistic timing parameters according to GR (e.g. §4.1 
of |Stairsl ( |2003] l). 

The to used in the analysis was the sum of the kinematic 
and relativistic terms held fixed while all other pulsar timing 
parameters (astrometric, rotational, and Keplerian orbital pa- 
rameters) were free to vary. We recorded the of each timing 
solution, assigning a probability to each grid point based on the 
difference between its and the global minimum. We used 
this ensemble of probabilities to calculate confidence regions 
for the parameters of interest and to place limits on the stellar 
masses. The results are given in Figs. |5]|6] and |7] 

Figure |5] shows the orientations of the orbit allowed by the 
timing solution. The figure is a projection of the 95.4% confi- 
dence volume of the three dimensional analysis grid onto the 
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two dimensional space shown. There are four regions allowed 
by the timing data. The Keplerian orbital parameters, combined 
with the lack of Shapiro delay restrict /, but because the Shapiro 
delay depends on sin /, the resulting constraint is degenerate: if 
i is allowed, so is 180° - /. Since there is a detectable i , Eq.[T] 
restricts 0^ - Q to two possibilities (one positive and one nega- 
tive value) for any given value of /. Thus there are four regions 
of allowed solutions, labeled A, B, C and D in the figure. 

Figure |6] shows 68.3%, 95.4%, and 99.7% confidence lim- 
its on cos / and Mj space. For reference, the figure shows lines 
of constant mass difference mc - nip, at intervals of Q.2Mq , as 
calculated from / and Mj, and it indicates the region at low 
inclination angle that is excluded because mp > is not sat- 
isfied in that region. The confidence limits were calculated by 
marginalizing over the probability values for values of 0^ - Q 
for a given combination of cos / and Mj. 

The best-fit solutions have relatively high - nip, i.e. high 
companion masses and low pulsar masses, although the confi- 
dence contours stretch to much lower values of mc - nip. The 
determination of Mi is dominated by the relativistic 6j , but the 
obtained values are perturbed by the kinematic Sio , given by 
Eq. [2] This perturbation splits the allowed values for the total 
mass into two regions. This can be understood by plugging the 
values of 0^ - f2 from each of the four regions of Fig. |5] into 
Eq. 2, and noting that esc i is always positive, so that dco is 
positive for solutions A and B and negative for solutions C and 
D. This means that the observed 66j has been biased towards 
higher values for solutions A and B, so that the true relativistic 
6<jj is lower than that calculated without the kinematic correc- 
tion. Since Mi is proportional to 6<jj^^^, this means that the 
true total mass is lower for solutions A and B. Similarly, the 
true mass is higher for solutions C and D. 

Figure[7]shows the 95.4% confidence volume projected into 
several two-dimensional parameter spaces. For this figure, the 
values of / and 0^ - Q correspond to solution D, but the dis- 
tributions of all other quantities are essentially identical for all 
four solution regions. We note that, for completeness, all pos- 
sible values for the parameters are shown in Figs. [5] |6]and[7] 
however we consider it unlikely that the pulsar mass will be 
lower than IMq . 

We used the probabilities from the grid analysis to con- 
strain nip and m^. The results are essentially identical for all 
four solution regions. The central 95.4% confidence intervals 
are nip = 0.72+°-5i and = 2.00+|J-^^ Mq . If, instead of 
central confidence intervals, we use 95.4% confidence upper 
and lower limits for the individual masses, these convert to 
trip < 1.17 Mq and ni^ > 1.55 Mq . Note that the apparent dis- 
crepancy between these numbers originates from diff'erent ar- 
eas covered in the probability distribution. For the total mass, 
the 95.4% confidence intervals are 2.7188 + 0.001 1 Mq for so- 
lutions A and B, and 2.7217 + 0.0018 Mq for solutions C and 
D. However, allowing any of the four solutions yields the range 
2.12Q^^^^^j^Mq . This last number is the most accurate value for 
Ml, as it reflects the uncertainty as to which orientation of the 
orbit is correct. 



4.2. Evolution 

The system is valuable for studies of DNS evolutionary sce- 
narios. It has been believed for a long time that supernova ex- 
plosions result in high kick velocities on the remaining neutron 
star ( ,Bailes||I989[ ). Recently, it has been argued that different 
evolutionary scenarios can be possible in binary systems with 



particular mass and chemical properties ( Podsiadlowski et al. 
2004 Van den Heuvel 2007 ). The so-called electron-capture 



collapse of an O-Ne-Mg core of a helium star leads to a/asf su- 
pernova explosion which is believed to be symmetric and there- 
fore does not result in a kick of the second-born neutron star. 
For most DNSs in which proper motion measurements have 
been possible, a low space velocity has been derived. Our mea- 
surement of the proper motion of PSR J1518H-4904, implying 
a velocity of about 25 kms"', is consistent with the evolution- 
ary scenario mentioned above, and another piece of evidence 
that not all pulsars receive a large kick at birth. 

Van den Heuvel (2007) argues that the low velocity of 
DNS systems may be correlated with the masses of the second- 
formed neutron stars in DNSs being somewhat lower than nor- 
mal: 1.25(6)Mo . Our 95.4% limij^of > 1.55Mq appears to 
contradict this correlation. 

In case of a symmetric supernova the mass lost during the 
explosion can be calculated by the expression e = AM^^/Mi 
( Bhattacharya & van den Heuvel|199T i, where e is the eccen- 
tricity and Mj the total mass after the supernova explosion. For 
PSR J15 18-1-4904 this results in AMsn = O.68M0 . If we use 
our 95.4% lower limit on the companion of 1.55 Mq this leads 
to a core-progenitor mass of 2.23Mo , which is considered to 



be in the lower range for helium cores ( Dewi & Pols 2003 1, 
and the progenitor star must have had a mass of about 9Mq . 
As these stars are beheved to produce degenerate O-Ne-Mg 
cores ( [Podsiadlowski et al.|2004i [Poelarends et al.|2008| l, this 
is another indication that the companion neutron star must have 
formed in an electron-capture core collapse. The somewhat 
higher companion mass can possibly be explained by assum- 
ing extra fall-back during the supernova event of ~ 0.2Mo (e.g. 
[Fryer & Kalogera[2001,,Fryer^2007) . To allow for this amount 
of fallback the supernova explosion has probably been very 
weak. An electron-capture core collapse is the most probable 
explanation for the parameters of this system. However, con- 
straining the masses even better would be very interesting. 

4.3. Search for the companion 



Apart from the original discovery papers [Nice et aT] ( |1996[ ); 
Sayer et al. ( 1997 1, there is no report of searches for pulsa- 



tions of the companion of PSR J1518H-4904 . The Green Bank 
Northern Sky Survey, optimized to find millisecond pulsars, 
had a flux-density limit at 370 MHz of 8 mJy for slow pulsars 
(with periods above 20 ms). 

The very high characteristic age and the low magnetic field 
imply that PSR J 15 18-1-4904 is the first-born and recycled neu- 
tron star in the system. The (unseen) second neutron star will 
be the young object, with a long spin period, which has proba- 

•* The 99.7% limit is m,. > 1.39Mo 
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bly already slowed down significantly and possibly passed the 
death-line to become undetectable as a radio pulsar Moreover, 
if the second neutron star is still active as a radio pulsar, it has 
a reasonable chance of being beamed away from us. In this 
case, taking into account that the precession timescale of the 
spin axis is very long, it is not expected that the companion, if 
active as a pulsar, will rotate into view on a short timescale. 

However to be certain that no weak pulsed signal from the 
second neutron star has been missed, we have searched several 
WSRT obsei-vations of 30 minutes at 840, 1380 and 2300 MHz 
for pulsations of the companion of PSR J1518-(-4904 . We per- 
formed an acceleration search on the data although, because of 
the wide orbit, the expected smearing due to the orbital motion 
of the companion is likely to be negligible. 

In the only double pulsar system known so far, 
PSR J0737-3039A/B, the second pulsar is only visible for a 
small part of each orbit ( jLyne et al.|2004| l. To take this possibil- 
ity into account, and assuming the second pulsar would be vis- 
ible in an equal fraction of the orbit as PSR J0737-3039B, we 
have searched another set of observations spread equally across 
the full orbital phase range. In all searches no evidence of pul- 
sations of a companion were found to the limits of 1 . 1, 0.25 and 
0.5 mJy at 840, 1380 and 2300 MHz respectively, which means 
that if the companion is active as a pulsar, and beamed towards 
us, it must be less luminous than 0.098 mJy kpc^ at 1380 MHz. 

Conclusions 

Using more than 12 years of timing observations obtained 
by 5 telescopes, we have improved the timing solution of 
PSR J15 18+4904. We accurately determined the total mass 
of the system to be Mj - 2.7183(7)M0 and the proper mo- 
tion fij - 8.55(4) mas yr"' . From our timing solution we were 
able to set constraints on the individual masses of the system: 
ffip = O.72+°-^iM0 and = 2.00+°-^^Mo (95.4%), and the in- 
clination of the orbit: / < 47 degrees (99%). 

Future prospects 

As discussed above, at the present timing precision it will not 
be possible to detect any further PK parameters in this system 
on any reasonable timescale. Recently we have begun observ- 
ing the system using the new PuMa II coherent dedispersion 
backend at the WSRT (Karuppus amy et al.|2008| l. These obser- 
vations use a total of 160 MHz of bandwidth and thus should 
have at least a factor V2 improved SNR. When combined with 
the somewhat narrower profile obtained using coherent dedis- 
persion (see FigjSjl, this should lead to improved arrival time 
measurements. At present we have insufficient data to check on 
the long term timing precision, however the error in the individ- 
ual arrival times are at least a factor two better than the previous 
data. These improvements might allow us to put stronger limits 
on the PK parameters in this system. 
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Fig. 2. Best fit timing residuals including all timing data avail- 
able. The top panel shows the different data sets, from bottom 
to top: GBT, Jodrell Bank, Westerbork, Effelsberg and Nangay. 
The offsets have been adapted manually to show the data sets 
separately. This plot was generated from the best timing solu- 
tion, which is shown in the middle panel as residuals vs time, 
and residuals vs orbital phase in the bottom panel. This tim- 
ing solution has an rms of 6 /^s and the parameters are listed in 
Table |2] 
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Fig. 3. Chisq contour plot of m^-sini parameter space. Lines 
with probabilities indicate regions with certain confidence lev- 
els. The dashed line represents the constraint on total mass 
resulting from the <1> measurement. From the non-detection of 
Shapiro delay we can set an upper limit to the inclination of the 
orbit: / < 47 degrees. This corresponds to a lower limit for the 
companion mass of 1 .29 Mq . 




Pulsar Mass (M^) 

Fig. 4. Mass-mass diagi-am for PSR J15 18-1-4904 and its com- 
panion. The hatched region is excluded for / > 90 degrees, the 
dash-dotted lines are indicating constraints resulting from the 
X measurement, / < 69 degrees and shapiro delay limit, / < 47 
degrees. The diagonal line constrains d>, with errors. The dot- 
ted lines for constant pulsar and companion mass indicate the 



range of neutron star masses measured in other DNSs (Thorsett 
& Chakrabarty|1999[m^ et al.|2004t|Faulkner et al.|2005| l. 
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Fig. 5. Allowed values for orbital orientation / and - D.. The 
figure shows a projection of the 95.4% confidence volume of 
the three dimensional grid analysis of timing solutions onto this 
two dimensional space. See the text for discussion. 
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Fig. 6. Allowed values for total mass and cos /. The two sets 
of contours are 68.3%, 95.4%, and 99.7% confidence limits on 
these quantities for solutions A and B (lower, white contours) 
and C and D (upper, gray contours). Dotted lines indicate val- 



ues of nir 



at intervals of 0.2Mq . The dark region on the 



right of the plot is excluded as it does not satisfy mp > 0. 
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Fig. 8. A pulse profile obtained using PuMa II at the WSRT. 
The data were obtained in a 30 minute observation at 1380 
MHz with a bandwidth of 160 MHz and were analysed using 
a coherent filterbank with 512 frequency channels and a time 
resolution of 6.4 /vs. A similar pulse phase range is shown as 
in Fig. [T] and the improved resolution shows somewhat more 
detail in the pulse profile. 



Fig. 7. Values of relativistic y, kinematic to, kinematic i, or- 
bital orientation 0^ - Q., and masses and nic allowed by the 
timing solution, as a function of inclination angle. The figure 
shows a projection of the 95.4% confidence volume of the three 
dimensional grid analysis of timing solutions onto each of the 
two dimensional spaces shown. The regions corresponds to so- 
lution D, but other solutions give essentially identical results. 



